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CO ' For the logarithmically singular parabolic equation (jl.ip below, we 

establish a Harnack type estimate in the L^^^ topology, and we show that 
the solutions are locally analytic in the space variables and differentiable in 
time. The main assumption is that In it possesses a sufficiently high degree 
. of integrability (see (|1.3p for a precise statement). These two properties 

are known for solutions of singular porous medium type equations (0 < 
m < 1), which formally approximate the logarithmically singular equation 
p.ip below. However, the corresponding estimates deteriorate as m — >■ 0. 
It is shown that these estimates become stable and carry to the limit 
as m — ^ 0, provided the indicated sufficiently high order of integrability 
is in force. The latter then appears as the discriminating assumption 
between solutions of parabolic equations with power-like singularities and 
logarithmic singularities to insure such solutions to be regular. 
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1 Main Results 



Wc continue here the investigations initiated in [2 |3], on the local behavior 
of non-negative solutions to logarithmically singular parabolic equations of the 
type 

u G aoc(0,T;LL(i^)), In " e LL(0, T; <f (i?)) ; 
ut - Ainu = weakly in Et ^ E x {0,T] 

where E is an open set in and T > 0. It is assumed throughout that 

u e L\^^{Et) for some r > max {l; ^} (1.2) 

and that 

Inu e Li^,(0,T;Lf^^(£;)) for some p > 1. (1.3) 

The modulus of ellipticity of the principal part is u^^. Therefore the equation 
is degenerate as u — > oo and singular as m — >■ 0. It was shown in [2] that 
p.2p implies that u is locally bounded in Ex, and hence the equation is not 
degenerate. Likewise if (|1.3|) holds for some p > N + 2, then the solution u is 
locally bounded below, and hence the equation is not singular. As a consequence 
u is locally, a classical solution to This was realized by establishing a local 

upper and lower bound on u, via a pointwise Harnack-type estimate. 

The main results of this note are that if m is a locally bounded, weak solutions 
to dTT]), then: 

i. If Inu satisfies (|1.3p for some p > 2, then u satisfies a local Harnack in- 
equality in the L^^^ topology, as opposed to a pointwise Harnack estimate. 

ii. If \nu satisfies (jl.3p for some p > + 2, then u is locally analytic in the 
space variables uniformly in t, and C°° in time. 



2 Harnack Type Estimates in the Topology of 

7-1 
^loc 

For /9 > let Kp be the cube centered at the origin of R''^ and edge p, and for 
y G let Kp{y) denote the homothetic cube centered at y. Moreover, Qp{9) 
denotes the parabolic cylinder Kp x (—0/9^, 0]. For < s < t < T and y ^ E let 
p be so small that K2p{y) x {s,t\ C Et- Since u G L^^{Et) the quantity 



M 



ess sup u 

K2pX{s,t] 



is well defined and finite. Also, if (|1.3p holds then the quantity 

u{x, t) 



Ap = ess sup 

S<T<t 



In- 



M 



dx 



(2.1) 



(2.2) 



is well defined and finite. 
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Proposition 2.1 Let u be a non-negative, local, weak solution to U.l\) satis- 
fying in addition lll.2\) and lll.3\) for some p > 2. There exists a positive con- 
stant 7 depending only on {N,r,p} and Ai and A2, such that for all cylinders 
K2p{y) X [s,t] C Et, there holds 

sup / ui.,r)d.<J inf / ui.,r)d. + l^\, (2.3) 

where 

\ = 2- N. (2.4) 



2.1 Weak Solutions Versus Distributional Solutions 

The L^^^ Harnack type estimate (|2.3p , continues to hold for merely distributional 
solutions to the second of (jl.ip whereby Inw is only in L\^^{Et)- The assumption 
(|1.2p however is still in force, and (|1.3|) is required to hold only for some p > 1. 
The constant 7 depends only on Ap for some p > 1. 

In §[5] we will prove (|2.3p first for such distributional solutions. The proof is 
rather simple due to the linearity of the principal part with respect to Inw. The 
linearity however is immaterial, as (|2.3p is a structural inequality valid for weak 
solutions to quasilinear parabolic equations with singularity and degeneracy of 
the same nature as (jl.ip . To be specific, consider non- negative, local, weak 
solutions to quasilinear parabolic equations of the type 

u&Ci,,{0,T;Ll,{E)), \nu e LI,{0,T;W,];^{E)) 
Ut ~ div A{x,t,u, Du) = weakly in Et- 

Here the function A : Et x R^+^ M.^ is only assumed to be measurable and 
subject to the structure conditions 

|p|2 

A{x,t,u,p) ■ p > Co 

" a.e. in Et, (2.6) 

\A{x,t,u,p)\<Ci^-^ 
u 

where Co and Ci are given positive constants. Assume that u and Inw satisfy 
pT ^ -pr ^ and introduce M and Ap as in (P7T|) and (j^ . Then u satisfies ([^ 
with 7 depending on the data {N, r,p, Co, Ci} and Ai and A2. 
The proof of this fact is more involved and it is given in § [B] 



3 Analyticity of Local Weak Solutions to (11.11) 

The precise statement of these results hinges on the notion of "intrinsic neigh- 
borhood" of a point (xo,to), as determined by the degeneracy and singularity 
of the equation in (|l.ip . 
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3.1 The Intrinsic Geometry of ( 11.11) and Main Results 
from [2] 

Let u be a non-negative, local, weak solution to Having fixed {xo,to) € 
Et, and Ksp{xo) C E, introduce the quantity 

0^=i^-(/ u'^{;to)dxy (3.1) 

where e S (0, 1) is to be chosen, and g > 1 is arbitrary. If 6* > assume that 

(x„, to) + QspiO) = KspM X (to - e{8pf,to] C Et. (3.2) 

These are backward, parabolic cylinders with "vertex" at {xo,to) whose height 
depends on the solution itself through the quantity 6. In this sense they are 
intrinsic to the solution itself. 

Continue to assume that u satisfies (|1.2p and let Inu satisfy (|1.3p for some 
p> N + 2. Then 

AI — ess sup u (3.3) 

{Xa,ta) + Q8p(e) 

is well defined and finite. Moreover the dimensionless quantity 

^JK.ix^)^ M ) \ KsMJ ^ ^ 

is well defined and strictly positive. Finally the dimensionless quantity 



Ap = ess sup 



is well defined and finite. 





u[x, t) 


p \ 




dx 1 




M 





forsomep>iV + 2 (3.5) 



Theorem 3.1 (Pointwise Harnack Estimate [2]) Let u be a non-negative, 
local, weak solution to satisfying the integrability conditions and 

il.Sp for some p > N + 2, and assume 6 > 0. There exist a constant e G (0, 1), 
077.6? a continuous, increasing function rj — > /(yy, Ap) defined m R"*" and vanishing 
atrj = Q, that can be quantitatively determined apriori only in terms of {N,p, q}, 
and Ap, such that 

inf m(-, i) > /(yy, Ap) sup u 

(3.6) 

for all i G (to - -^Op'^ , to] 

For 77 — >■ and Ap — > 00, the function rj — > f{rj,Ap) can be taken to be of the 
form 

f{i-l) = exp I - I /or < ?7 < 1 and Ap > 1 (3.7) 
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for positive constants Ci and C2 that can be determined apriori only in terms 
of {N, p, q} . Moreover 

£ — > and Ci + C2 ^ 00 as p ^ 00 or p ^ N + 2. (3.8) 

Remark 3.1 In [2 the constant i] was given a more general form. For the pur- 
pose of this note the definition (j3.4[) represents the degeneracy of the equation, 
quantified by M — > 00. The occurrence Ap — > cxd quantifies, roughly speaking, 
the singularity of the equation. 



3.2 Analyticity in the Space Variable, of Solutions to (II. ip 
at 

Theorem 3.2 Let u he a non-negative, local, weak solution to satisfying 
the integrability conditions lll.2\) and \1.3\) for some p > N + 2, and assume 
9 > 0. There exist two parameters C and H , that have a polynomial dependence 
on f{r]), [/(?/)]~^, N , such that for every N -dimensional multi-index a 

\D''u{xo, to)\ < . ' u{xo, to). (3.9) 

Moreover, for every non-negative integer k 

, N Ci/2fe(2fc)! , ,1 ,, ,^ 
-Q^u{xo.to) < -^-^uixo^of-''. (3.10) 

Remark 3.2 The theorem continues to hold, with the same assumptions, for 
local, weak solutions to the quasilinear equations (j2.5|) . provided the function 
A is analytic in all its arguments whenever u is bounded above and below by 
positive constants. 



4 Approximating (11.11) by Porous Medium Type 
Equations 

Consider local, non-negative, weak solutions in Et to the porous medium equa- 
tion 

ziGCioc(0,r;iL(^)), w^Lll^.T-Wl^liE))- 
ut- Aw = weakly in Et ^ E x {0,T] 

where 

u™ — 1 

w = for < 7Ti < 1. (4.2) 

m 

As TO -> 0, formally (|4.ip - (|4.2p tend to (|l.ip . In [3] a precise topology was 
introduced by which such a formal limit is rigorous. A natural question is 
whether solutions to (|4.ip - (|4.2p satisfy a version of the Ll^^ Harnack estimate 
(|2.3p . which as TO — >■ tends, in some appropriate sense to be made precise, 
to that of Proposition 12.11 A similar issue arises for the local analyticity of 
Theorem O 
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4.1 Harnack Type Estimates in the Topology of L^^^, for 
Weak Solutions to m]) - d4:2|) 

A first statement in this direction is that u satisfies 



sup / u{x, T)dx < 7 

s<T<tJKp{y) 



inf / u(x,T)dx 1 — r- 



(4.3) 



where 

A A^(m - 1) + 2. (4.4) 

Here 7 depends upon TV and m and 7(771) — > cx) as to — > 0. Thus, one cannot 
formally recover (|2.3p by letting to — > in (|4.3p . However, (|4.3p is rather 
general as it continues to hold for non-negative, local weak solutions to general 
quasi- linear version of (|4.ip . Precisely 



u e Cioc(0,r;LL(^^)), ^ e iL(o,T;<i(i?)) 



(4.5) 



ut — div A(a:;, t, u, Du) = weakly in Et 

where the function A : Et x R^+^ — > is only assumed to be measurable 
and subject to the structure conditions 



A(x,t,u,p) •p> Cou'"-'b|" 
\A{x,t,u,p)\ < Ciu"'-'\p\ 



, a.e. in Et, (4.6) 



where Co and Ci are given positive constants. In such a case the constant 7 in 
(|4.3p depends also on these structural constants. The proof of these statements 
is in [4], Appendix B. 

A major difference between (|2.3p and (|4.3|) is that in the latter u is not 
required to be locally bounded, nor does 7 depend on some analogue of the 
quantity Ap as defined in (|2.2|1 . This raises the question as to whether (|4.3p 
holds with 7 independent of m but dependent on some analogue of Ap. 

Henceforth we assume 

u e L[o(.(£'t) for some r > maxjl; ^(1 - ™)} (4.7) 

and that 

weL^,{0,T;Ll^{E)) for some p > 1. (4.8) 

It was shown in [3] that (|4.7p implies that u £ L^^{Et) and hence the corre- 
sponding quantity M defined as in p.ip is well defined and finite. Set 



A„,p = esssup / J ^ ] dx\ . 4.9 

This is the analogue of (|2.2p . and, if (|4.8p holds, it is well defined and finite. 
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Proposition 4.1 Let u be a non-negative, local, weak solution to the singular 
equations ^.5^ ~ (4^ j *^ -^t- There exists a positive constant 7 depending only 
on the data N^Co,Ci, and Ani^i, Air 2 o.i^'d independent ofm, such that ^.3^ - 
\4-4\ l holds true, for all cylinders K2p{y) x [s,t] C Eq^. 

As a consequence, (|2.3p can be recovered from (|4.3p . with the indicated de- 
pendences, as TO — > 0, provided proper conditions are placed, that insure the 
pointwise convergence of the solutions to (|4.ip - (|4.2p to solutions of These 
conditions are identified in and we will touch on them briefly in the next 
subsection. 



4.2 Analyticity in the Space Variable, of Solutions to (14. ip 



at (Xojto) 



Having fixed (xojto) G Et and K^p[xo) C E, the intrinsic geometry of (|4.1 
(|4.2p is determined by 



e[f u'i{x,to)dx] . (4.10) 



The intrinsic cylinders are as in p.2p with replaced by 6'^. The analogues of 
T] in (|3.4p are 

,u{x,to)Y 



where r > 1 is any number such that 



M 



dx 



i _2_ 

(4.11) 



A,. = iV(TO - 1) + 2r > 0. (4.12) 

In [4] a Harnack estimate of the form of (|3.6p was proved for these solutions 
with /(•) depending only on a and of the form 

/(^) = ^ (4.13) 
7(to) 

where 7(m) ^ 00 as m — t- 0. The constant /3 depends on and /3(Ar) — >■ 00 
as Ar 0. It was observed in [J § 21.5.3 that, for each fixed to S (0, 1), such 
an estimate implies the local analyticity of the solutions in the space variables 
about {xojto), and at least the Lipschitz continuity in time. 

Because of the indicated dependence of 7 (to) on to in (j4.13p . such a regularity 
does not directly carry to the limit as to — 0. In [3] we established a Harnack 
estimate of the form p.6p for solutions to (|4.iP " (|4.2p and its quasi- linear versions 
l|4.5p - (|4.6p . with /(•) depending on cr, as defined in (|4.1ip . and A^ p as defined 
in (|4.9p provided p> N + 2. The form of such /(•) is the same as that in p.7p 
with the proper change in symbolism. The new feature of such an /(•) is that, 
while depending on the quantities a and Am,p, each quantifying the degeneracy 
and the singularity of the equation, is independent of m and hence is "stable" 
as TO — >■ 0, provided a and A^.p arc uniformly estimated with respect to to. 
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As a consequence, the analyticity estimates of Theorem l3.2[ can be recovered 
from the analogous ones for solutions to (|4.1[) - (|4.2|) whenever solutions {um} 
to the latter converge pointwise to solutions to the former. In [3] it was shown 
that this occurs if there exists m,, £ (0, 1) such that 



u„ 

Wr. 



e Li^^(0, T; L\^SE)) for some r > max{l; ^N} 
, e (0, T; Ll^{E)) for some p>N + 2 



(4.14) 



uniformly in m £ (0,m**). It is also required that there exists an open set 
Eo C E and a positive number <JEo;T such that 

/ Um{-,T)dx > (Je„:T uniformly in to. (4-15) 

5 Proof of Proposition [2JJ for Distributional So- 
lutions to (11.11) 

The proof is a local version of an argument of [1] for global solutions to the 
porous medium equation for < ?7i < 1. Let Ci £ C^(M^) be such that 

< Ci < 1 

Ci = 1 in Kp (5.1) 
Ci = in R^\K2p. 

Then, by the divergence theorem 

ACidx= [ ACidx= [ ^ds = 0. (5.2) 

By (|5.2p . for any positive constant M, any as in (j5.ip . and any non- negative 
function v such that A^i Inw is integrable, we have 

/ ACilnuda;= / ACiln( — )da;. (5.3) 
Now consider C,2 G +oo) and C,i as in (|5.ip . such that 

p p 

By the previous notation, with u a solution to (jl.ip . we have 

/•OO /• /"OO /• 

— / / (^'2(^iudxdt — / / ^2^Ci Inudxdt. 
Taking into account (|5.3p . for any positive constant M 



Ciudx^ [ ACiln(^)dx in l?'(0,r), 
dt Jk2„ JK2„ ^1 
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and also in Lj^^{0,T). Therefore 



dt 



^ / Ciudx = f ACiln(^)dx < / lACil In(^) 



dx. 



By the definition of Ap and from the previous estimate, 



dt 



C^iudx 



K2 



< Ar. 



\ACi\Pdx] \K, 



K2 



. C(N) jv N C(A„,N) 

p^ p'^ 

where A = 2 — N. Taking into account the size of the support of Ci, for any 
0<s<i<rwc conclude 



u{x,t)dx < C{Ap, N) i / u{x,s)dx-\ — 



6 Proof of Proposition 12.11 for Weak Solutions 

to (^:e)-(^M) 

6.1 An Auxiliary Lemma 

Lemma 6.1 Let u be a non-negative, local, weak solution to the quasi-linear 
singular equations \2.5\ )- ^2^) . in Et ■ There exist two positive constants 71, 72 
depending only on the data {N ,Co,Ci} , such that for all cylinders K/^pijj) x 
[s,t] C Et, and all a G (0, 1), 



^C'rf^d^<7i(l + Ai)5, + ^(A? + A2) ' 



where 



Sa = sup / u{-,T)dx. 

s<T<tJK(i + „^piy) 



Proof - Assume {y,s) = (0,0), fix cr g (0,1), and let x — )■ ({x) be a non- 
negative piecewise smooth cutoff function in K(^i^„-^p that vanishes outside 
i<'(i+o-)p, equals one on Kp, and such that \DC,\ < (ap)^^. Let si £ [0,t] be 
such that 



Sa = sup / u{-,s)dx^ I 



w(-, si)dx. 



We also set 



^ def 
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In the weak formulation of (|2.5p - (|2.6p . take the test function 



and integrate over Q = K(^i^^^p x (0,i], to obtain 
0=JJ ^u(^n—^ C^dxdT + J J A{x,T,u,Du) ■ D 

= h+l2- 

We estimate these two terms separately. 



In 



So 



dx dr 



h 



d 



In 



So 



C^rfa; dr 



S, 



Qn[M<5„ 



dr 



S. 



u In — C dx dr 



A'(i+„)pn[«<S„] 



-f (i+<T)pn[u<5„ 



C\x) [u\n — + u]{x,t)dx 
C{x) I u\n—+u] {x,0)dx 



Next, 



In 



A(a;, r, u, Du) ■ D 

A(a;, r, u, Du) ■ D 

Qn[ti<5„] 



dx dr 



u 



dx dr 



Du 



C A{x,T,u,Du) dxdr 

Qn[u<Sa] ^ 

C In — A(x, r, u, Du)DC,dx dr 

C3n[«<5„] 



< -Co 
+ 2Ci 



- ^ 2 

Qnlu<s„] " 



dx dr 



^l^^\M\DC\dxdr 



Qn[ti<5<,] 



Qn[«<5„] " 



u u 
e\^dxdr^ ^ 



2 „2 



// - 




/ JQn[u<s^] 


u 



where 7 = Therefore, we conclude that 



Co 

2 



/ ^ 2 

iQn[u<S„] u 



2 ' ^ dxdr < 



-f (i+<T)pn[M<5„] 



dx dr. 



C'ix) [ u In — + u ] {x, t) dx 
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7 



// - 




/ jQn[u<s^] 


u 



dx dr 



= J1 + J2. 
We have 

Ji - 



< 



-f (i+<T)pn[M<5„] 
f 

-f (i+<T)pn[M<5„] 



C (x) I u In h u I (x, t) dx 



u In — + u j (a;, t) dx 



S.I ^iJ-^Mdx^l uMdx 

< Sa [ In — (x, t) (ix + Scr 

JK(i+„)pn[«<5„] " 
<lSA \iY—{x,t)dx + S. <7Ai5„+5^ =7(1 + Ai)5^, 



where 7 = 2^. Moreover, 



cr^P^ JjQn[u<s^] 



< 



7 



u 



dx dr 



dx dr < 



-2 „2 







// - 




/ jQn[u<S^] 


u 



dx dr 







f'l 




Jo JK,, + „y^ 


u 



dx dr 



< 



sup 







/ 






u 



Hence, we have 



iQn[u<s„ 
Now, if we take the test function 



C'l^dxdr < 7(Ai + 1)5. + ^Al(-^ 



(6.1) 



in the weak formulation of (|2.5P " (|2.6p and integrate over Q = A'(i_|_cr)p x {0,t], 
we obtain 



^u^hi^^ (^'^dxdT+ II A{x,T,u,Du) ■ D 
h + h- 



dx dr 
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We estimate these two terms separately. 



Next, 



Qn[M>5„ 



-f (i+ff)pn[«><s„] 



C {x) yu. In J uj (x, dx 

(x) In J u'j (x, 0) dx 



A{x, T, u, Du) ■ D 



In. 



Qr\[u>Sa] 



A{x, T, u, Du) ■ D 

Du 



dx dr 



dx dr 



C A{x,T,u, Du) dx dr 

Qn[«>5<,] " 



+ 2 

>a 



u 



In ^ A(x, T, u, Du)DC,dx dr 



Qr\[u>Sa] 

2 dxdT 

Qn['u>5„] ""^ 



> 



2Ci 



Qn[M>5„ 



Qr\[ii>Sa 



Cl^JL\R^\DC\dxdr 

] Oct U 

e'^^dxdr ^ 



2 „2 







// - 




1 jQr\[u>s^] 





dx dr. 



where again j ~ 2-^. Therefore, we conclude that 



Co 

2 



< 



C 2 

Qn[ii>5„] " 



^ ^ ^ da; dr < 



if(i+„)pn[u>5„] 

1 



C i) dx + 



-'^(i+<T)pn[ti>s„] 
7 



win--- {x,Q)dx 



t2 o2 







// - 


In^ 


1 JQr\[u>S„] 


5. 



dx C?T 







// - 




/jQn[M>5^] 





< 



-R'(i+,T)pn[u>5„] 
< + ( In 1^ I 5, 



C{x) uln-^ {x,Qi)dx + S, 



dx dr 



C{x) ( uln^ ) {x,Q)dx + S„ + 



In- 



M 



tp 



N 



In- 



M 
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We need to evaluate \nM/Sa- As in the interval (0, 1] the function /(s) = — Ins 
is convex, Jensen's inequality yields 

ln^=ln. ^-m/ 



Sa ii.^ u(x,Si)dx Jt^, , M 



, , u(x,Si) , f , M , 

< i - In ^ ' dx ^ f In — rdx 



< 7A1, 

where 7 = 2^^. Hence, we have 

C^lE^dx dr < 7(Ai + 1)5. + ^ A? f 4r V (6-2) 

The lemma follows by combining estimates (|6.ip and (|6.2p . 

The use of ^In as test function can be justified using ^In and 

then letting e — > 0. ■ 

Corollary 6.1 Let u he a non-negative, local, weak solution to the singular 
equations i2.5\) - l[276\) . in Et- There exists a positive constant 7 depending only 
on the data {N,CotCi}, such that for all cylinders Kip{y) x [s,t] C Et, and all 
ae (0,1), 

1 



\A{x, T, M, Du)\dx dr 
<7max{(l + Ai)^;(A2+A2)3| 5^ + J_(^^) 



t - s 
7^ 



Proof - Assume (y, s) = (0,0), and let Q ~ Kp x (0,t]. By the structure 
conditions of A 



— / / \A{x,T,u,Du)\dxdT < — 2- (( J ^-dxdr 

pJo Jk„ P JJq u 



By Lemma 16.1 1 we conclude. 



6.2 Proof of Proposition 12.11 

Assume {y, s) = (0, 0). For n = 0, 1, 2 . . . set 

E" 1 ~ Pn + Pn+1 J-, 
— p, Kn = Ap„; p„ = , Kn = ilp 

j=l ^ 
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and let x Cn{x) be a non-negative, piecewise smooth cutoff function in Kn 
that equals one on Kn, and such that l-DCnl ^ 2"+^ /p. In the weak formulation 
of (|2.5p - (|2.6p take as a test function, to obtain 

I u{x,Ti)(„dx < / u{x , T2)(ndx -\ / / \A{x,T,u,Du)\dxdT 



< I «(a:,r2)C„dx + 7 2"5„%i - +7 4" - , 



where 



<S„ = sup / u{-,T)dx. 

0<T<t J Kn 



Since the time levels ri and T2 are arbitrary, choose T2 one for which 



u{-,T2)dx~ inf / u{-,T)dx 



K2n 



def 



I. 



With this notation, the previous inequality takes the form 

t 



5„ <Z + 7(data,Ai,A2)4"(4) + 7(data, Ai, A2) 2"5i+i(4^) 



By Young's inequality, for all So G (0, 1) 

<Sn < £oSn+i +7(data,Ai,A2,£o)4'' 

From this, by iteration 

5o<e;'5„+7(data,Ai,A2,eo)[l+ (^]] "e(4£o)'. 

Choose £0 so that the last term is majorized by a convergent series, and let 
n — 00. ■ 
The proof of Proposition 14. 1 1 for weak solutions to the porous medium type 
quasilincar equations ()4.5p - (j4.6p . is similar, with the obvious modifications, and 
we confine it to Appendix B. 



7 Analyticity in the Space Variables, of Solu- 
tions to (11.11) 

Let M be a non-negative, local, weak solution to (jl.ip . satisfying the integrability 
conditions (|1.2|) and (|1.3|) for some p > N + 2. 

Fix (xo,to) G Et, assume that Ksp{xo) C E, and assume that the quantity 
6 defined in (|3.ip is positive. The cylinder {xo,to) + Q%p{d) is assumed to be 
contained in the domain of definition of u as in (|3.2p . The quantities M, r] and 
Ap are defined as in p.3p - p.5p . 
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From the Harnack-type inequality p.6p . 

Ifiv)] u{xo, to) < u{x, t) < [f{ri)]~^u{xo, to) 
for any (x, t) within the cyhnder 



(7.1) 



Q = K2pixo) X (to- —ep^,to]. 



Equation (|l.ip can be rewritten as 



(7.2) 



ut - diY{ - Du) =0. 



(7.3) 



By (j7.ip this can be regarded as a particular instance of a hnear parabohc equa- 
tion with bounded and measurable coefficients. By known results (for example, 
0, Chapter II) local, weak solutions to (|7.3p are locally bounded and locally 
Holder continuous. Consequently, (|7.3p can be regarded as a linear parabolic 
equation with bounded, and Holder continuous coefficients. Again by classical 
theory (see [7], Chapter HI), one can conclude that local, weak solutions are 
indeed C°° with respect to the space variable. 
By (|7.ip the quantity 9 can be estimated as 



9<e sup u{-,to) < e[f{ri)] ^u{xo,to) 
e>e inf u{-,to) > e[f{?])]u{xo,to). 



def 



Let S ~ S{ii) = e[f{ri)] and introduce the change of variables 



t-to 



It maps Q onto to 



A def „ , 1 
Q ^ K2X { 



u[xo,to)p^' u{xo,to)' 



ir ,^,0]DQ/=i^A-,x(-A,0], 
lbu{Xo,to) Id 



and within Qs the function v satisfies 



1 \Dv\^ 
vt = —, 



(7.4) 



(7.5) 



(7.6) 



with 

f{v)<v<f{ri)-\ 
By a result of [6], there exist constants < cr < 1, C and H such that 

suplZJ^vl < Ci?l"l|a|!, 
< Ci72fe(2fc)! 



sup 



(7.7) 
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where Qas = x {—jqO-^iO]. Tracing the dependence of constants gives 

C = 7iCo, H^^2^MCo[.f{v)]-\[f{vT^} (7.8) 

where 71 and 72 are constants independent of v and Co is a function of f{r/) 
and satisfies 



Qk 

dt' 



N 

<Co inQ,- for |a|+2fc<4[y] 



16, 



where [a] denotes the integer part of a. Thus in particular an upper bound on 
these derivatives up to the indicated order, gives their analyticity as signified by 
(|7.7p . Assuming such an upper bound for the moment, we return to the original 
coordinates to get 



\D''u{xo,to)\ = |Z?"i'(0,0)| 



u{Xo,to) 



Qk 



«(0,0) 



(7.9) 



„2k 



„2k 



The proof is concluded, once the dependence of Co on /(yy) is determined. This 
estimation can be achieved by local DeGiorgi's or Moser's estimates. While the 
method is known, it is technically involved and reported in detail in Appendix A. 

The analogous analyticity estimates for solutions to the porous medium type 
equation (|4.ip - (|4.2p are similar, with the obvious changes, and we omit the 
details. 



Appendices 

A Analyticity in the Space Variables, of Solu- 
tions to ( fTTI ). Estimating the first 4[f ] + 16 
Derivatives of v. 

Wc will use expressions such as w'^f(w), f and similar ones, but wc only have 
at our disposal the notion of weak solution, and therefore, such a way of working 
docs not seem justified. However, by the Harnack estimate of Theorem 13. H 
solutions arc classical, and in these calculations wc arc turning the qualitative 
information of u being classical into the quantitative information of u being 
analytic. 

With respect to the previous sections, we use a different notation for cylin- 
ders, and we let Q{p, 6) = Kp x (—0, 0]. 
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A.l An estimate of ||Dt;||oo 

Take D^. of the logarithmic diffusion equation and set Wi = v^.. to get 



dtWi — div i-Dwi -WiDv) = 

V V 

Setting w = {wi, . . . , wn), yields 

Wt — div {-Dw \rw ® w) =0. (A.l) 

For all derivations below we stipulate that A<l,A = A^^>l,A<i<A, 
6* < 1 and p <\. 

Proposition A.l Let w be a solution to II A . 1]) and be a cutoff function in 
Q^Q{p,9). Then 

sup / / sf{s)ds('^ dx 
-e<t<oJKp Jo 

+ ^ JJ^\Dw\^fi\w\)e dxdt+^ JJ^\w\\D\w\\^n\w\)C'dxdt 
^^11 fi\w\)\w\'\DC\'dxdt + ^ ff \w\'fi\w\)edxdt 



A JJq a jjq 

^JJq 



\wff'{\w\)C^dxdt+ // / sf{s)ds2CC,tdxdt, 
JJq Jo 



where f : R+ — > R+ is a bounded, non-negative, non- decreasing Lipschitz con- 
tinuous function. 

Proof - Multiply (jA.ip by wf{\w\)C'^ and integrate by parts to get 



= J J — J sf{s)dsC'^dxdt + I 
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where 
1 = 



= JJ (^^Dw - ^w®w^ ■ [Dwf{\w\)C +wf'{\w\)D\w\(^ 
+wf{\w\)2CDC] dxdt 

-f(\w\)\Dw\^C^ + -wf(Uv\)Dw ■ D\w\(^ + -wf(\w\)2CDw ■ D( 

V V 

Dwf{\w\)C^ - —w (g) w ■ wf'{\w\)D\w\C^ 



;W ®w ■ wf{\w\)2(DC 



dxdt 



> 



-\Dwff{\w\)e + -\w\\D\w\fr{\w\)e + -\w\f{\w\)2CD\w\ ■ DC 

V V V 



^\w\wD\w\f{\w\)e - ^\w\'wD\w\f{\w\)e - ^\w\^wD(:f{\w\)2C 



Observing that iDjiyH < \Dw\ the previous inequahty yields 



dxdt 



< 



- sf{s)dse + MDM'f{\w\)e + XMD\w\\'n\w\)e 
[A\w\f{\w\)2CD\w\ ■ DC + A^\wfD\w\f{\w\)C^ 



dxdt 



A'\w\'D\w\f{\w\)C^ + A^\w\^\DC\fi\w\)2C] dxdt 



< 



A 
2 

A4 
A 
A 
2 



\Dw\^f{\w\)C^ dxdt + 



A2 



f{\w\)\w\^\DC\^dxdt 



w\'^f{\w\)C^dxdt + A'^ // \wf\DC\f{\w\)2Cdxdt 



w\\D\w\\^f'{\w\)C' dxdt + 



Q 



Q 
2A 



wff{\w\)C' dxdt. 



Wc win use this energy estimate and Moscr's iteration to derive a bound for 
\\Dv\\oo. 

Proposition A. 2 There exists a positive parameter 5 that depends only on 
such that if p & {^i^] then 



where 7, pi, and fi2 positive parameters that depend only on N. 
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Proof - We first estimate ||w||2- In the weak formulation of (|1.1[) take the test 
fmiction w^^. By standard calculations 

A// \Dv\^C^dxdt<j I v^C^dx + jjl v\DC\^ dxdt. 

JJq{p,b) JKf,x{-e} JJQ{p,e) 



Hence 



\Dv\^ dxdt < 7 

Q(ap,ae) 



A2 



A (l-a)2 



1 + - 

p . 



i < 7 T 



1+1 



XJ (l-a)2p2' 

(A.2) 

The energy estimate (jA.ip with /(|w|) = |i(;|2^ and /? > 0, yields 



\w\^^+Wx,t) dx + ff \Dw\''\w\^^C''dxdt 
^^11 \w?^^^m\' + Q] dxdt + ^ [[ \w\^^+'e dxdt. 



Q(pfi) ^ •l-lQi.pfi) 

Notice that v is locally Holder continuous with the Holder norm and its exponent 

' A 
A 

|«(x,-)~«(0,-)| <7Ap° 



a depending only on y and N , that is 



Now we apply an integration by parts and use the Holder continuity of v to 
estimate the last term of (|A.3[) . If we freeze the time variable, then 



[ \w\^^+^C^dx^ [ \w\^P+^wwedx 
|^|2/3+2^2^^, • D{v - v{0, •)) dx 

{v - v{0, •)) [Av\w\^^+^C^ + 2\wf^+^CDC ■ Dv 
+ {2(3 + 2)\w\'^^+^D\w\DvC'^] dx 
< 7A / [Av\w\^^+\^ + 2\w\^^+\DC ■ Dv 

+ {213 + 2)\w\^^+'^D\w\ ■ DvC"^] dx 

<l ( |«;|2^+4(2rfa; + 7p'"A2 / \w\^P\Dw\^dx 

4 JKp JKp 

+ 7p2"A2 / \w\^^+^\D(:\^ dx. 
Jk„ 



(A.4) 
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Since all constants arc independent of t, combining (jA.3[) - (|A.4p yields 

sup / \w\^^+^C^{x,t)dx+ ff \Dw\^\w\^^C^dxdt 
-e<t<oJK„ JJQ(n.e) 



(A.5) 



A2 . 

+ 1^ II P^" [\w\^^\Dw\^ + \w\'^+^\DC\'] dxdt. 
Thus, by taking 

we have 

sup / \w\'^^+'^C^{x,t)dx+ II \D'w\'^\w\^^(^ dxdt 

-9<t<oJKp JjQipM) 

II \w\'p+'[\Dc? + Q+p''^m^]dxdt 

II \w\'^+^[\DC?+Q]dxdt, 
as p^" < 1 by (lA.SP . An application of the Sobolev parabolic embeddings gives 



// 



9) I A JJq{pM) 



Now take 

Pn^(jp+ ^-^^P, 9n=Cr9 + ^-^0, Qn = Q{Pn, On), 

let ^ be a standard cutoff function in Qm and set 

/iV + 2\" ff 
a„ = 2(^ — — — j and ^ J J \M°'"dxdt. 

Let 6 = 4° and a = 1 + Begin from 2(3 + 2 ~ Uq, that is from /? = 0, and 
apply the above estimate recursively up to 2/3 + 2 — a„ . This gives 

W<C°6"/,:, where ^^-^{(r^+(T^ 
Iterating these recursive inequalities yields 

In+l < C^"=o''"^'6^-«^"-*)"7f 
< (^7(7V)a" + i^7(Ar)a" + ija" + i^ 
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Now take the power of both sides and let n — oo. Taking into account the 
estimate of Iq in (|A.2p . we obtain 



\\w\\oc,Q{ap,a0) 



(i-ay 



lo 



where the constants /ii and fi2 have been properly modified. 



A. 2 An Upper Bound for ||t't||oo 

Proposition A. 3 Letv be a classical solution to the logarithmic dijfusion equa- 
tion and assume < A < v^^ < A in Q{p,9); then 

/Ay^(i + 0^ 

where 7, fix, and fj.2 are positive parameters that depend only on N. 

Proof - Multiply p.ip by the test function wjC^ and integrate over the cylin- 
der Qi — ^+°'^ ), where a G (0, 1). Here ( vanishes on the parabolic 
boundary of Qi and takes value 1 in Q{ap,ad). A standard calculation gives 







vfO^dxdt + -Dv[C^Dvt + 2vtCDC] 

V 



dxdt 



v1('dxdt+^('^\Dv\^ + -VtCDvDC 
2v at V 



dxdt 



Vf(^'^dxdt 



\Dv\ 



This gives the estimate 



K p + crp X{0} 
1 , 1 

-^vt^ + -CCt 

2v V 



1 

2^ 



iDvl'^C^ dx 



dxdt 



-VtCDvDCdxdt. 

V 



II v'i(^dxdt< II \Dv\'^^\vt\(^dxdt+ 1 1 -C\Ct\\Dv\'^ dxdt 
+ 11 -\Dv\\vt\2C\DC\dxdt 



1 

< - 

- 2 



v'^C'^dxdt + 



Qi 



\Dv\ 
4^4 



■Cdxdt 



'i\Dv\ 



\DCrdxdt 



\Dv\' 



■\Q\dxdt. 
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Taking into account the estimate for Iji^w ||oo,Qi of the previous section, we have 



/Ay^(l+^ 

\vth.Q(.p.a0)<l[jj TT^^ 



for some /ii(A^), and H2{N) > 0. 

Now take the time derivative of the logarithmic diffusion equation, and in 
the corresponding weak formulation use the test function vtf{\vt\)C'^ where / : 
K+ — s> R+ is a bounded, non-decreasing Lipschitz function, and C vanishes on 
the parabolic boundary of Q = Q{p,9) and takes value 1 in Q(ap, a0). Let 
M = \\Dv\\ oo,Q- A standard calculation yields 

JJq l§i\''t\^f'^^dxdt + ^JJ^ \Dvt\^jedxdt 

+ xjj \Dvt\^\vt\f'C^dxdt 

<2A jj^\Dvt\\vt\K\DC\dxdt 

+ A^JJ^ \vt\\Dv\ [\Dvt\fC^ + \vt\f'\Dvt\e + '2\vt\fC\DC\] dxdt 
= Ii+ I2 + h + h. 
Let us estimate the four terms. 

h<\jj^ \Dvt\^fCdxdt + ^ JJ^ \vt\y\DC\^dxdt; 

h<^jj^ \DvtffCdxdt + 4^— j \vtffC^dxdt; 

/3< 2 jj^\Dvtf\vt\f'edxdt + 4^— jj^\vtff'Cdxdt; 

l4<MA^JJ \vt\^fC^dxdt + MA^ JJ \vt\^f\D(:\^dxdt. 

Summarizing we have 

sup / / sfis)dsC^dx + X \Dvt\^fC^dxdt 
-e<t<oJKf, Jo JJq 

+ X JJ^\D\vtW'f'\vt\Cdxdt 

< J^)\m^ + 1] ff [\v,\'fe + \vt\'re + H'fiDcn dxdt 



2 I j j sf{s)dsC(tdxdt. 
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Now take /(s) = sl^ for /3 > 0; then 

sup / \vtf+^C^dx+ ff \D\vt\\^\vtfC^dxdt 

P + Z ~9<t<aJKf, JjQ 

< l(^jJ[M^ + 1] JJ^ + \DCf + 101] dxdt. 



Let a = l + jf, 



A 



(l-CT)2p2 {l-a)e 



def 



and let w = Vt- An application of the Sobolev embedding yields 



\wm^)^dxdt 



Q(ap,ae) 

< ( sup / \\'w\^C\'^dx 
KCil + Pf^-l^jj \wf+^dxdt^ . 



\D[\w\ — (:]?dxdt 



Take Q„ as before and define 

+ 2 f N + 2' 

/3o = 0, /3„+i+2 = (/3„ + 2)^^ ^ 2(^—1 -2, 



N 



N 



and 



It then follows that 



\P"+^dxdt. 



In+l < C'h^I^ 

for some positive constant h depending only on N . A standard iteration gives 

^ ^7(Ar)a" + i^7(Ar)Q" + ija" + i^ 

Therefore, taking the — ;rfT power of both sides and letting n — > oo, we have 

\\vt\\oo,Qi.p,.0) <C'^b'' JJ^\vt\'dxdt. 

To conclude this section, bound the right hand side using the estimate for ||wt||2 
on an intermediate cylinder. ■ 
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A. 3 An Upper Bound for ■^D°'v 

Differentiating the logarithmic diffusion equation successively we have 



Qk 



1 „ d'' 



]<k 



jj df^-i V dV 



DD°'v 



V dp 



D^Dv 



(A.6) 



0, 



|/3|<|a| ^'"^^ j<k 

where fc G N and a is a multi-index. For an integer n > let 

Qk 



E 

k-\-\a\—n 



at'' 



We have the following bound for general derivatives of the logarithmic diffusion 
equation. 

Proposition A. 4 Let v be a classical solution to the logarithmic diffusion equa- 
tion in Q(p^6) and fix a G (0,1). Assume < A < v^^ < A in Q(p,9). There 
exists a positive parameter 6 that depends only such that if p E (0,(5], then 
in Q{a-p, <j9) 



(1 +6)-M2) 



(A.7) 



(1 - (jY^ 

where 7, pi and p2 are positive parameters that depend only on N and n. 

Proof - Muhiply (TOI) by the test function ^D°'vf{\w\)C where / : R+ 
IR+ is a bounded, non decreasing Lipschitz function. Here C vanishes on the 
parabolic boundary of Q2 = Q{ (^+'^)^ ) and takes value 1 in Q{ap,a9). 

Standard calculations and a sum over fc + |q;| = n give 

d /■'"'' 



Q2 



dt 



I sf{s)dsdxdt + \ // \w\\D\w\\^f{\w\)C,'^dxdt 
J JQ2 



A E 

A;+ |a|— n 



\D-^D''v\^f{\w\)edxdt<I, 



where 

/ = - 



E 

A;-f |a|— n 



lu{\M)^D-vD^D-vDCdxdt 



- E 

fc+|at = 
l/3|<l"l 3 

d 



E 

3<k 
k 



jJ dt^-1 V dP 



DD°'v 



J J dt''-^ V dP 



d' 



Qk 



I{\M)C;D^D"v + C^ — D"vf'{\w\)D\w\ + 2U{\w\)^D^vD(: 



dt 



dt 



dxdt. 
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Notice that 



E 

j<k 



ij df'-J V dV 



\l3\<\a\ ^' ' j<k 
<-P[|w| + l] 



,1 



for some polynomial P with variables {||^£'"z;||oo,Q2 for fc+ \a\ < n; j}. Thus 
E // \D^,D-vff{\M)edxdt 



k--\-\a\—7i ' 



A2 
2A 



E 



fe+|Q| 



D"v\'f{\w\)\DC\^dxdt 



1 II \D^D-v\'f{\w\)edxdt 



fe+|a| 



P^{\w\^ + l)f{\w\)C'^dxdt 



\D\w\\^f'{\w\)\w\Cdxdt + ^ // P^(\w\^ + l)\w\f'{\w\)Cdxdt 



P{\w\ + l)\w\f{\w\)C\DC\dxdt. 

Q2 

Wc obtain the following energy estimate 

.\w\ 



sup 

2 

A 



.tl^<t<0JK^ JO 



sf{s)dsC'^dx+^ [[ \w\\D\w\\^f'{\w\)C^dxdt 
2 JJoo. 



E 



<7 



A 
A 



\D^D'^v\'fi\w\)edxdt 



P^{\w\^ + l)fi\w\)Cdxdt+ // |w|V(h|)|£'CI 



Q2 



+ 2 



s/(s)(isC|CtMa^c?i, 



Q2 "'0 
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and also 



sup 

A 

^2 



:f{s)dsedx+^ JJ \w\\D\w\\'r{\w\)C 



dxdt 



\D\w\\^f{\w\)C'^dxdt 



<7 



P^{\wf + l)f{\w\)C^dxdt + 



\w\^f{\w\)\DC\^dxdt 



+ / / P^{\w\^ + l)\w\f'{\w\)C^dxdt 



sf{s)ds (^\Qt\dxdt. 



Q-i Jo 



Now we have at our disposal the sup-estimates for Dv and Vt in terms of 
A/A only. Next we assume the supremum of -^D^v is estimated for all 
fc + |a| < n by a similar quantity as the right hand side of (|A.7|I , By Moser's 
method, the above energy estimate will yield a bound for the case fc + |a| = n. 
These will depend on the norms of ^D°'v and a polynomial with variables 
{||^D«7;||oo,Q. forfc+|«| <n;A}. 

Take f{s) = for /? > 0; then the energy estimate yields 



1 



sup 



/3 + 2 _e+a^^^^QJK^ 

3 



\wf+\''dx + 



\D\w\\'^\wfC'^dxdt 



<7 



1 



1 



1 



1 



(l-CT)2p2 {l-a)e 

Let a = 1 + f and C = {^f 
then an application of the Sobolcv embedding yields 



1/3+2 



\w\'^]dxdt. 



and assume \Q2\ < 1; 



// 

J Jc 



Q{ap,ae) 

< I sup 



||W| 2 (l^'dx 



1 /3 + 2 



+ l)dxdt 

\w\^+^dxdt) +C"p2«(l+^)3a^ 



after a proper adjustment of the constant 7 in the definition of C. Take 



Pn = 0-/9 + 



1 - cr 

2n+l ' 



2"+i 



Define 



/3o = 0, + 2 = + 2) 



iV + 2 
N 



N + 2 
N 
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and 

/r, = / / \w\^"+^dxdt. 



Wc have that 

In+i < C"P^°6"/° + 

where the constant b depends only on N. A standard iteration and a proper 
adjustment of P give 

In+i < 2^"=i(""~i'(cp2)^"=i'"'6^"=o("-')"7;j"''' 

Take the — ttft power of both sides and let n — > oo to obtain 

lklloo,Q(.p..e) < {CP'r f f \w\^dxdt + {CP^y. (A.8) 

Remember that Q2 — s±a9_y^ order to conclude the proof of Propo- 

sition IX4I we only need to estimate ||w||2,Q2- 

A. 3.1 An Estimate of ||w||2,g2 

It is enough to give an estimate of ||w;||2, Q(o-p,cr6»)- Replace k in (|A.6p by — 1 
and assume fc > 1. We can rewrite (|A.6p as 

Qk-l Qk-l \ 

5^^"--div -i.^Z.".).div/ 



where 



If C is a smooth function in Q2 and takes value 1 in Q{crp, ad), then the standard 
estimate for the linear parabolic equations gives 



\-g^D"^'\\2,Q{c,p,ae} 



< 



1 A..^'-..,„. ^ 



Let us denote P as a polynomial of variables {||^-D"w||oo,Q2 f^i' < n-] j}- 

Observe that 

\\divf\\lQ^<jP[h+h + l3 + l] 
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where 



h 



h = 



d 



k-l 



Qk-2 



Dv 



dxdt: 



gk-l 



D 



dxdt: 



dxdt. 



Here for an integer I 



'„,|2 



\0\=i 



These quantities can all be estimated in the same way. Indeed, they all contain 
spatial derivatives and we can use the principal part of the differentiated (|f .ip 
to estimate them. Precisely, we write ()A.6P as 



c\k-s / 1 fl*^-^ 

— D^vt - div ( -D^:^^D'^v + 

j<k-s 



dt*" 



E 

\f}\<\ri\ 



E 

j<k-s 



V dt'' 



J J dt^"-"-'-' V dv 



0, 



DD'^v 



where 1 < s < fc and \ri\ = \a\ + s. Take the test fimction 

Qk — s 



where C, vanishes on the parabohc boundary of Qa = 3(i+£)0^ 
takes value 1 in Q2- Integrating in Qa, a standard calculation yields 

-\D—^D'^v\^C^ dxdt 



d_ 

, dt 

E 

j<k-s 



:ik—s 



-D'^v 



dxdt 



dt''-' 
k - s\ d''-"'^ 1 



j J dt''-"-^ V dti 

E f;0 E rDSS--^^-'- 



l/3|<|r,| 

Qk 



j<k-s 



D 



dtk-' ^ ^ dtk 



j J dt''-"-^ 

Qk^S 



:D''vCDC 



< 



\D^^^D^v\^edxdt + ^ 



p. 
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This together with (jA.9[) gives 

ll^^"«ll2,Q(ap,.e) <7 



1 



1 



{l-a)e (l-a)V 



P. 



(A.IO) 



Now the only remaining case is \a\ ~ n. For this, we consider the equation 
(|A.6p with fc = and assume \a\ = n — 1. Take the test function C,^D°'v, where 
C vanishes on the parabohc boundary of Q3 and takes 1 in Qiap, a9). 



-\DD°'vYCdxdt = 

V 



-DD"vD°'v2CDCdxdt 



jja~f} _DD^\>DD°'vC'^dxdt 



|/31<|a|"'"'^3 

V // (Z]D''-^-DD'^vDD°'vD"2CDCdxdt 



<^ [[ \DD°'v\^(^dxdt^ 
Summing over all jal = n — 1 actually gives 



{l-(T)e (1-ct)2/,2 



P 



E 

|a| — n 



\D°'v\'^C'^dxdt < 



1 



1 



(1 - cr)2p2 



P. 



If we take into consideration an intermediate cylinder, then this, together with 
((XTO| in (|X8l) . yields 



kllc 



•.Q{cyp,cre) 



< pi 



1 



(A.ll) 



for some 7 depending only on N . The induction hypothesis and the definition 
of p in (jA.5|) imply that 



/A V 1 + e^^^ 



(A.12) 



B Proof of Proposition 14.11 for Weak Solutions 
to Equations (l43|)-(l4:6|) 



An Auxiliary Lemma 

Lemma B.l Let u be a non-negative, local, weak solution to the singular equa- 
tions OT^-OTfil), in Et- There exist two positive constants 71, 72 depending 
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only on the data {N,Co,Ci}, such that for all cylinders K^p^y) x [s,t] C Et, 
and all a G (0, 1), 



\^(^^dxdT < 7i(l + A^a)p'^^ Sl-' 



2^2 V"^^,! 



where 



Scr = sup / u{-,T)dx. 



Af(l- 



Proof - In the following we restrict to < m < | , since we arc mainly interested 
in proving the stability of the estimates as to — > 0"''. For to G (j, 1) similar 
arguments hold, provided a slightly different test function (p is chosen (see [3], 
§ B.1.1 for more details). 

Assume (y, s) ~ (0,0), fix <t G (0,1), and let x — )• ({x) be a non-negative 
piecewise smooth cutoff function in K(^i_^_„-'jp that vanishes outside /•C(i_|_ct)p, 
equals one on Kp, and such that \D(\ < {ap)~^. Let si G [0,t] be such that 



Scr = sup 

0<s<t 



/ 



m(-, s)dx 



u(-, Si)o?x, 



-R'(i+<T)p(i/) 



and set 



^ def Sa 



In the weak formulation of (|4.5p - (j4.6p take the test function 



and integrate over Q = K/i_f_^\p x (0, t], to obtain 



d I u 
97" - 



A(x, T, u, Du) ■ D 



C dxdr 



U 2 -Sa 



dx dr 



h+h 
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We estimate these two terms separately. 

,/,/o 9t" 1 m 



C dx dr 



d I u '2 ~ So- 
u 



Qn[u<5„ 



dr 



u{x,t) 



(^dx dr 



-s; 



■ ds \ dx 



ds dx. 



Next, 



I2 = A{x,T,u,Du) ■ D 



Qn[u<s^ 



A(a;, T, u, Du) ■ D 



/ + 
m 



dx dr 



dx dr 



< 



< 



1 

'2 
2 

a 
2 

2Ci 

' 4 

7 



Qn[u<s„ 



Qn[u<5„ 



Qn[u<5^ 



C^w 2 ^ A{x,T,u, Du) ■ Dudx dr 



\ m 

u-'^-^u"'-^C^\Du\^ dx dr 



A(x, T, w, Du) ■ DC, dx dr 



u 2 ^ Sa 



u"'-^\Du\\DC\dxdT 



Qn[«<5„] \ 

2„.-S--2| 7-,„.|2 



Qn[M<5„] 



Cu-^-'lDwI^da; dr 



<^^P^ J jQn[ii<s„] 



2i 

M2' 



U 2 -Sa 



dx dr. 



where 7 = Therefore, we conclude that 



Ca 

Co 

4 
< 



Qn[u<5„ 



Cu--^\Du\'' dx dr 



7 



U2' 



(2;,0) 



■ ds \ dx 



dx dr 



m 
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= Jl + J2. 



Wc have 

Ji - 



s 2-5^ 



■ ds \ dx 



I 



< I s; 

/if(i+„)pn[«<5„] 



^ ( ^ I \dx 



< 



< 



r 

-f<'(i+^)pn[M<5„ 
f 

-Ff(i+<T)pn[M<5„ 
f 

f 

-R'(i+<T)pn[u<5p 
1 



5^ 



y ^ -^ y 

ym 2 — yjiy ™ 



5^ 



2 - m 
2 

2 - m 



dx 



w(x, 0) 



1-^ 



- 1 



i(x,0) 



2 it(x,0) / fu{x,0) 



2 — m 2 — m Sa \ \ S, 



- 1 



1 2 u{x,0) 



2 — m 2 — m iS^ 



2 - m 



5o- ^ dx 



-'<'(i+,T)pn[M<5„] 



^ u(x,0)— ^— ^ — ^ — dx = J{ + j". 



2 — m 



-'<'(i+<T)pn[«<5„ 



m 



J[ = ^—Sl " / dx<^— 5;^ "p^^, where 7 

2 - m ^K(i+„,pn[u<5.] 2 - m 



2 — m 



2 — m 



-R'(i+<T)pn[u<5„] 



-»(x,0) 2 2 
m(x, L)j dx 



-i_^u(x,0)i-^5.=^ -m(x,0)^ , 

Oct ; — — dx 



-f<'(i+<T)pn[«<5p] 



5^ 



< ^5^-- 
2 — m 



- u(x, 0)' 



-dx 



'if(i+„)pn[«<5.] mM2 

<-^5r~A.^.i/o'^ = :7^Azf,i5^"^p^^, where 7 = 2 

' 2 — m 
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2 — m 
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Therefore, 



Ji <7(1 + A™,i)p^^5^ ' 



Moreover, 

2 

dx dr 



2r 
U2 



- / / " ^ ' "~ — - 1 dx dr 



< -^Sj I I ( 1 dxdr 





5j 


- 




mM 









2 



2 

<^Y^Sjtp^ sup + ( , ,1* ) 



Cr2p2 2 



Hence, we have 

Co 



4 JjQn[u<So 



C'^u'^-'^\Du\^dxdT 



Now, if we take the test function 



? 9 



2 



(B.l) 



in the weak formulation of (|4.5p ~ (|4.6p and integrate over Q = K^i^^^p x (0,t], 
we obtain 



A(a;, r, u, Du) ■ D 



'Q 



+ 



dx dr 



33 



We estimate these two terms separately. 



Qn[u>s„ 



u — 



dr 



Next, 



I4 ~ A{x, T, u, Du) ■ D 



(^^dx dr 



m 



ds dx 



ds dx. 



Qn[M>5„ 



A(.T, T, u, Du) ■ D 



/ + 
Sg -u 
m 



dx dr 



dx dr 



> 



> 



+ 2 

Co 

2 

-2Ci 

Co 



Qn[u>5„ 



C^u 2 ^A{x,T,u,Du)-DudxdT 



Sa ^ - U 



Qn[u>s^] 



A(x, T, u, Du) ■ DC,dx dr 



Qn[u>5„ 



Cu-^-'^u'"'-'\Du\''dxdT 



u'^-'^\Du\\DC\dxdT 



Qn[u>s^ 



Qn[u>s^ 



Cu--^\Du\^dx dr 



7 



o'^P^ J jQn[u>s„ 



U2' 



dx dr. 



p2 

where again 7 = 4-^. Therefore, we conclude that 



Co 

4 

< 



Qn[u>s„ 



Q^u-~^\DuYdxdT 



-R'(i+<T)pn[u>5„ 
7 



.(.,0) §- 



■ ds\ dx 



U2' 



^ P JJQr\[u>s„] \mu^Sa 



dx dr = J3 + J4,. 
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We have 



i(x,0) 



■ ds ] dx 



< 



-f (i+a)pn[«>5„] 



u{x,0} 1 



dx 



< 



< 



f 

-R'(i+CT)pn[u>5„ 

-R'(l + <r)pn[u>5„ 

r 

-R'(i+CT)pn[M>5CT 
1 



1 -i i-i 



TO TO 



TOl/' 



2m 
? 

2 - to' 



dx 



u{x,0) 

5a 



m(x, 0) 



2 — TO \ y iSct 



2 u(x,0) / /m(x,0) 



2 — m 2 — m iS„ 



2 w(x,0) 



2 - m 
2 



2 — TO 2 — TO 5a 



iSct ^ dx 



5a 

u{x,0) 



TO 



-R'(i4.„)pn[«>5„] 



2 - TO 



uix,0 f^' -"("'°)"" dx = 4 + 4- 



2 - TO 



dx<—^sl "p^^ where 

2 — 771 



2 - TO 



5^ = - u(x,0) 2 
u(x, 0) dx 



< 



Mi + ")p 
2 5a ^-M 



-5. 



2 — TO TO 2 — TO toM 2 

As in the interval (0, A/] the function 



2 M^-5a^ iv^ci-^ 
P 2 5a 



/(«) 



mM' 



is convex, we can apply Jensen's inequality and conclude that 



2 

2 — to' 



J3 < 



J K, 



M 2 — u(a;, si)" 

TOMt 



-dx 
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Therefore, 

J3 <7(1 + A^,i)p'^^'5^^. 

As for J4 we have 



7 1 / - 5^ 



2 



< -Tr-^-=— — t sup / u^™'(x,t) dx 



3, 



2 



^2^2 5m I 



^0<T<tJK 



Co 

4 



where we have taken into account Jensen's inequahty once more. Hence, we 
have 

r r 

C^w^-2|Dupda; dr 
'Qn[u>5„] (-g 2) 

< 7(A^,i + + 4^A|. ,Sjtp^('-^\ 

2 (T^p^ 2 ' 

The lemma follows by combining estimates (jB.ip and (|B.2[) . 

The use of (^ " ~:^f'^ ^ ^ as test function can be justified using u + e 

instead of u, and then letting e — > 0. ■ 

Corollary B.l Let u be a non-negative, local, weak solution to the singular 
equations ( |^.5[ )-( f77^ , in Et ■ There exists a positive constant 7 depending only 
on the data {N, Co-, Ci}, such that for all cylinders K4p{y) x [s, t] C Et, and all 
ae (0,1), 

I f ,., , 7, .9 .9 , 1 / 1 — s 



Aix,T,u,Du)\dxdr<^A%^,+A%,,)^^ S:^ 

1 

+ 7(1 + A^.i)^(^j S.^ 
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Proof - Assume (y, s) = (0,0), and let Q = Kp x (0,i]. By the structure 
conditions of A 

- [ [ \Aix, T, u, Du)\dx dr <— [[ u"'^^\Du\dxdT 
PJo Jk, P JJq 



< 



Ci 



<7 



Du I dx dr 



(1 + A^,iy^^5, 



u^'^dxdr 



9 ? 



t sup 

0<r<tJK^i+„^p 



u2™(a;, T)dx 



<7 



(1 + A^,i)p^^5^ ' 



By simple computations, we conclude. 



Proof of Proposition 14.11 

We conclude as in the proof of Proposition l2.H relying on Corollarv lB.li instead 
of Corollary 16.11 ■ 



References 

[1] M.A. Herrero and M. Pierre, The Cauchy problem ut = Au™ when < 
m < 1, Trans. Amer. Math. Soc, 291(1), (1985), 145-158. 

[2] E. DiBenedetto, U. Gianazza and N. Liao, On the Local Behavior of Non- 
Negative Solutions to a Logarithmically Singular Equation, Discrete Con- 
tinuous Dynamical Systems Ser. B, 17(6), 2012, 1841-1858. 

[3] E. DiBenedetto, U. Gianazza and N. Liao, Logarithmically Singular 
Parabolic Equations as Limits of the Porous Medium Equation, Nonlinear 
Analysis Series A: Theory, Methods & Applications, 75(12), 2012, 4513- 
4533. 

[4] E. DiBenedetto, U. Gianazza and V. Vespri, Harnack's Inequality for De- 
generate and Singular Parabolic Equations, Springer Monographs in Math- 
ematics, Springer- Verlag, New York, 2012. 

[5] E. DiBenedetto, Y.C. Kwong and V. Vespri, Local Space-analyticity of 
Solutions of Certain Singular Parabolic Equations, Indiana Univ. Math. 
J., 40(2), (1991), 741-765. 



37 



[6] D. Kinderlehrer and L. Nirenberg, Analyticity at the boundary of solutions 
of nonlinear second order parabolic equations, Comm. Pure and Appl. math 
XXXI, (1978), 283-338. 

[7] O.A. Ladyzenskaya, N.A. Solonnikov and N.N. Ural'tzeva, Linear and 
Quasilinear Equations of Parabolic Type, Translations of Mathematical 
Monographs, 23, American Mathematical Society, Providence, RI, 1967. 



38 



